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We analyze the process of simultaneous photon pair generation and quantum walks realized by 
spontaneous parametric down conversion of a pump beam in a quadratic nonlinear waveguide array. 
We demonstrate that this flexible platform allows for creating quantum states with different spatial 
correlations. In particular, we predict that the output photon correlations can be switched from 
photon bunching to antibunching controlled entirely classically by varying the temperature of the 
array or the spatial profile of the pump beam. 

PACS numbers: 42.50.-p, 42.82.Et, 42.65.Lm 



Spontaneous parametric down-conversion (SPDC) is 
probably the most commonly used process for generation 
of quantum correlated photons pL] with many applications 
including quantum cryptography Q and quantum logic 
devices [3, 4]. However, the use of bulk optics for gener- 
ating correlated photons as well as for the building blocks 
of logic gates hinders the scalability of the quantum cir- 
cuitry with increasing number of components. Indeed, 
the successful operation of a quantum optical circuit re- 
quires that the fidelity of the quantum interference, which 
lies at the heart of single-photon interactions, is preserved 
after passing through all optical components. Integrated 
optical quantum circuits have been seen as the solution 
for on-chip scalable quantum networks with important 
demonstrations of multi-photon entanglement [5], quan- 
tum factoring algorithms [6], and polarization entangle- 
ment [7|]. Additionally, integrated circuits are compact, 
stable, and could lead in the near future to mass produc- 
tion of chips for quantum computation. 

A particularly important platform for quantum ma- 
nipulation and integrated control of non-classical light is 
provided by an array of coupled optical waveguides [8]. 
Waveguide arrays have been used to perform quantum 
walks of photon pairs resulting in nontrivial quantum 
correlations at the array output [9, 10]. Such kind of cor- 
related walks involving quantum interference of several 
walkers can provide a speed-up of quantum algorithms 
delivering an exponential acceleration with the number 
of correlated walkers 12 |- However, in all schemes to 
date, the correlated photon pairs were generated exter- 
nally to the array by using bulk photonic elements. Such 
bulk elements may introduce quantum decoherence and 
impose stringent requirements on the losses associated 
with the connection of the array to the photon sources. 

In this Letter, we propose and demonstrate numeri- 
cally a novel scheme of quantum walks, involving simul- 
taneous generation of correlated photon pairs through 
SPDC and their quantum walks inside a single photonic 
element - an array of quadratic nonlinear waveguides. 
This scheme avoids entirely the need for complex inter- 



faces required in previous experiments [lo| , but most im- 
portantly enables novel ways for control of the spatial 
quantum correlations at the array output. In particu- 
lar, we show that by simply varying the phase-matching 
conditions for the SPDC process or the spatial profile 
of the pump beam it is possible to control the output 
quantum states incorporating photon bunching or anti- 
bunching. Importantly, such unusual quantum statistics 
are not possible when the photon pairs are created ex- 
ternally to the array. Although integrated photonic cou- 
plers [13, 14] and circuits [15] incorporating SPDC were 
proposed previously, we emphasize that integrated SPDC 
and quantum walks in a single nonlinear array lead to 
additional quantum interference between probabilities to 
generate photon pairs in different places of the array. 

Arrays of quadratic nonlinear waveguides are widely 
explored for manipulation of optical pulses at telecom 
wavelengths through cascaded generation of second har- 
monic [8]. Here, we consider the reverse process of SPDC 
and study the generation of correlated photon pairs as 
illustrated schematically in Fig. [IJa). To demonstrate 
the flexibility in controlling photon states, we consider 
a phase-matched near-degenerate type I SPDC, when a 
pump beam generates photons of the same polarization 
and frequencies approximately half of the pump beam fre- 
quency. Whereas non-degenerate SPDC can also occur in 
the array, it can be excluded through frequency filtering 
at the output. In the array, the signal and idler photon 
states can be associated with^extended Bloch waves of 
the form exp(z/c^^n + iPs,iz) [8], where n is the waveg- 
uide number, k;^- and Ps,i are the normalized transverse 
wavenumbers and propagation constants for signal and 
idler waves, denoted by s and i subscripts, respectively. 
The quantum walks in an array can occur due to photon 
tunneling between waveguides, and we consider a com- 
mon case when such tunneling occurs between neighbor- 
ing waveguides and its rate can be characterized by the 
coupling coefficients Cs,i @, Eoj. We assume here that 
the coupling coefficients are the same for signal and idler 
photons within the filtered frequency range, and denote 
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FIG. 1: (color online) (a) Schematic illustration of a quadratic 
waveguide array: the pump beam generates photon pairs 
that couple to the other waveguides, (b) Propagation con- 
stant vs. the normalized transverse wavenumber for the near- 
degenerate signal (red dashed line) and idler (blue solid line) . 
(c) Phase mismatch vs. frequency detuning for degenerate 
type I oo-e SPDC in a single LiNbOs waveguide of L = 5 cm 
for extraordinary narrow-band pump at 582 nm and ordinary 
signal and idler waves. 



Here 5 is a constant, is a pump normalized 

transverse wavenumber, n is a waveguide number, 
(cjmin, ^max) IS the filtered wavelength range, a de- 
fines pump spectrum, a) are photon creation opera- 
tors at the specified transverse wavenumbers and fre- 
quencies, |0,0) is a vacuum state. The function Vt de- 
fines an overlap between Bloch waves, Vt (/c^, /c^, k^) = 

En {kp) exp [i{k^ -k^- ^^z^)^] , where Ak{k^) is a 
/c-space pump spectrum, found as a Fourier-transform of 
the input pump profile An{n). The function ^ in Eq. (j2j) 
defines the phase-mismatch as ^{k^ ,k:^ ^Ug^uJi) = 
sinc(A/3L/2) exp(— 2A/3L/2), where the first term defines 
phase-matching width, the second terms accounts for 
the phase accumulated by traveling photons fv?\ and 
A(3{uJs,uJi) = p^^\ujs^uji)-p{ujs,k^)-p{uji,kl), where 
P^^^ {ujs^ooi) is a propagation constant of the pump beam. 

To be specific, below we consider a continuous wave 
or narrow-band pump with central frequency uo^^ , which 
spectrum can be taken as a(ijJs +cji) = 8{uo^^ —ujs—uji). 
Then, the final expression for two photon state is: 



C = Cs^i' Then, the propagation constants follow a gen- 
eral relation for waveguide arrays [8] , 

= P{uj,,ukti) = + 2Ccos(fci;,). (1) 

Here (3^^^ is a propagation constant for a single waveg- 
uide, and ujs.i are frequencies of signal and idler waves. 
We plot the characteristic dispersion curves for the signal 
and idler waves in Fig. [IJb). The dependence of prop- 
agation constant on transverse wavenumber gives rise 
to quantum walks of photons [lO|, which rate is higher 
for larger C. For a pump beam with optical frequency 
ujp '2uJs,i, dispersion relation analogous to Eq. ^ will 
also apply, however the corresponding coupling coeffi- 
cient Cp would generally have a much smaller value com- 
pared to the signal and idler waves, Cp <C C, due to the 
weaker mode overlap between neighboring waveguides at 
higher frequencies [8[. In practice, CpL <C 1, where L 
is the array length and therefore the coupling effects can 
be neglected for the pump beam by setting Cp = 0. In 
this case, the input pump beam profile An{n) remains 
constant inside the array. 

Now we study photon correlations at the array output 
by adopting the mathematical approach of Ref. [l6|, . 
We describe the photon states at the array output that 
are associated with the extended Bloch waves, which 
complete set can be defined by the transverse wavenum- 
bers from the first Brillouin zone, — tt < A:^ < TT. Then, 
the two-photon quantum state can be written as: 

xn{k^, kf, ki) $ [ki , fc,^ , w7t^,) 
xa{ujs +u}i)a^ (ujs,kj-) [uji.k^) |0,0). 



dk^dk^- / dAw\^k{k^,kl,ALu)}, 

(3) 

where 

\^k{k^,k^, Alo)) = Ak{k^ + ki) sinc(A/3L/2) 
eyi-p{-iAI5L/2)a){Aoj,ki^)a){-Auj,k^)\Q,Q). ^ ' 

We define the phase mismatch using Eq. (|T|), 
Ap{ki,kl,Auj) = AP^°\Alo) - 2Ccos{k^) - 
2Ccos 

{k:t)^ where Af3^^\Auj) is a mismatch in a sin- 
gle waveguide, and Auj = uJq — uj^^ /2 = uj^^ /2 — cj^ is a 
signal/idler frequency detuning from the degenerate fre- 
quency uj^^ /2. 

The phase-matching for degenerate type I SPDC can 
be achieved using temperature tuning in birefringent 
crystals with second order nonlinearity, such as LiNbOs. 
Waveguide arrays in LiNbOs can be routinely fabricated 
by Ti-indiffusion [1, IH, EH, whereas to reduce photore- 
fraction from the strong pump beam [19i] it is possi- 
ble to use MgO-doped composition [205 high sam- 
ple temperatures. Characteristic dependance of phase 
mismatch on the frequency detuning Ap^^\Auj) for the 
type I oo-e SPDC in near-degenerate regime is shown 
in Fig. mc). This dependance can be approximated by 
the second order polynomial, Ap^^\Auj) = Ac})^^^ / L — 
A(j)^'^'> Auj'^ /[{ujf^YL]. Here A^^^^ is a constant phase 
mismatch and Acj)^^^ corresponds to quadratic dispersion. 

We now calculate the second-order correlation function 
Tk{ks^ki) = J dAcj(|(^/c|^/c)p), which defines correla- 
tions between photons with specific transverse wavenum- 
bers. In order to determine correlations in real space 
with waveguide numbers and rii for signal and idler 
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FIG. 2: (color online) Correlations of photon pairs (a,c,e) in k 
space (spatial spectrum) and (b,d,f) in real space (waveguide 
numbers) for a pump coupled to a single waveguide n = 
with phase mismatch (a,b) A^^^^ = 0, (c,d) Acj)^^^ = 40. 



FIG. 3: (color online) Correlations of photon pairs (a,c) in 
/c-space and (b,d) in real space for A^*^^'* = and a pump 
coupled to waveguides n = 0, 1 with amplitudes (a,b) An{0) — 
A^{1) = 1 and (c,d) An(0) = -A^{1) = 1. 



photons, we apply Fourier- transform to the two pho- 
ton state |^/e(A;^,fc,^, Acj)) |^n(ns, n^, Acj)). We 
then calculate the photon number correlation function 
Tn{Tis^ni) = JdAcj(|(^n|^n)P)7 which can be mea- 
sured by two detectors scanning across the waveguide 
array output and measuring the coincidences [10| [see 
Fig. [Ha)]. We take parameters for a LiNbOs waveguide 
array [18] LC = 10 such that the total structure contains 
about 6 coupling lengths, and consider frequency filter- 
ing Acj^^nA^(i)/4°^ = -2 to Aa;^axA0(i)/4°^ = 2, 
which corresponds to several nm wavelength filtering for 
a 5 cm long array and a pump wavelength 582 nm. Then 
the required integrating time window for the detectors 
should be of the order of several hundred fs or longer. 

We first consider the case when the pump is cou- 
pled only to the central waveguide n = 0, which cor- 
responds to a constant Fourier spectrum of the pump, 
Ak{k^ -\- k:^) = 1. In Fig. [2] we plot photon correlations 
at the array output in /c-space and real space for different 
values of the phase mismatch A^^^^ in the range from 
to 40 (corresponding to heating the LiNbOs array from 
200 to 202. 5°C). For A^^^) = we observe a square shape 
for the /c-space correlations [Fig. [2fa)]. This indicates a 
pronounced correlation between the generated signal and 
idler photons with transverse wavenumbers satisfying the 
relations kj- ± k^ o:^ ±7r, which appears because these 
wavenumbers correspond to most efficient phase-matched 
interactions with A/3 = —2C cos [k^) — 2C cos [k:^) = 
at Auj = 0. We note that according to the disper- 
sion relations ([1]), the propagation direction can be cal- 
culated as [8] u{k) = -dpL/dk = 2CLsm{k). Then, 
we see that at phase-matching the photons in a pair 



would have the same or opposite propagation directions, 
i^{k^) ±:v{k^). Indeed, the corresponding real-space 
correlations shown in Fig.[2fb) reveal that the probability 
of detecting signal and idler photons in either the same 
waveguide (n^ = n^, bunching) or opposite waveguides 
{ris = — n^, anti-bunching) is significantly higher com- 
pared to the other probabilities. When A^*^^^ is increased 
to 40, only the generation of photons with kj- ~ /c^^ ~ 
takes place [Fig. Efc)], which destroys signal and idler 
correlations in real space [Fig. [2fd)]. 

The output photon statistics can be tailored by chang- 
ing the pump profile and phase. When the pump beam is 
coupled with equal amplitudes and phases to two neigh- 
boring waveguides, An{n) = 1 for n = 0, 1, then the k- 
space correlation pattern is strongly modified compared 
to pump input to a single waveguide [Fig. [3fa)]. This 
happens because the pump spectrum is primarily con- 
centrated in the central part of the Brillouin zone with 
\kp\ < 7r/2, and since for phase-matched interactions 
kp k^ -\- k^^ this suppresses generation of photons 
with wavenumbers k^ + k^ ±7r. The remaining phase- 
matched processes with kj- — kf- zLtt correspond to 
opposite velocities of generated photons since i^{k^) ^ 
— v(k^)^ and accordingly the anti-bunching regime pre- 
vails in real space [Fig. [3fb)]. When we introduce a tt 
phase difference between pump amplitudes in the input 
waveguides, i.e. ^^(0) = 1 and ^^(0) = exp(z7r) = — 1, 
then the situation is effectively reversed with the other 
phase-matched processes dominant in /c-space, leading 
to bunching of signal and idler photons in real space 
[Fig.lHcd)]. 

The output photon statistics can also be controlled 
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FIG. 4: (color online) Correlations of photon pairs (a,c) in 
A;-space (spatial spectrum) and (b,d) in real space (waveguide 
numbers) for a pump coupled to a two waveguides n = 0, 2 
(a,b) and n = 0, 3 (c,d). 



by couphng the pump to spatially separated waveguides. 
The correlation distributions look especially interesting 
for the cases when waveguides n = 0,2orn = 0,3 are 
pumped. The output /c-space interference patterns and 
real space bunching or antibunching in this case depends 
on whether the number of non-pumped waveguides in 
between the pump inputs is even or odd, c.f. Fig. HJ^a- 
d). This is analogous to the behavior previously observed 
in waveguide arrays with photon pairs coupled from ex- 
ternal source [ol, fioj, however in the case of nonlinear 
waveguide arrays with combined SPDC and quantum 
walks, the correlations seem much more significantly pro- 
nounced due to the quantum interference between prob- 
abilities to create photon pairs in different places along 
the pumped waveguides. 

In conclusion, we have studied simultaneous SPDC and 
quantum walks in an array of quadratic nonlinear waveg- 
uides and shown that the output correlations can be ef- 
fectively controlled by changing the relative phase of the 
pump in the two input waveguides, as well as by altering 
the phase mismatch for the degenerate type I SPDC pro- 
cess. Such control can enable careful engineering of the 
output quantum state, including dynamic switching from 
anti- to bunching regimes. Our results provide new av- 
enues for the development of quantum integrated circuits, 
combining generation of photon pairs and simultaneous 
transformation of the correlated photon states. It may 
be of interest to analyze similar systems with modulated 
profiles of the waveguide coupling tailored to specific re- 
quirements, or to consider type-II SPDC and generation 
of four-photon states ^] . 
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